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Abstract. Deep cxtragalactic surveys with ISOCAM re- 
vealed the presence of a large density of faint mid-infrared 
(MIR) sources. We have computed the 15 /im integrated 
galaxy light produced by these galaxies above a sensitivity 
limit of 50 /iJy. It sets a lower limit to the 15 /im extra- 
galactic background light of (2.4 ± 0.5) nW m -2 Hz -1 . 

The redshift distribution of the ISOCAM galaxies is 
inferred from the spectroscopically complete sample of 
galaxies in the Hubble Deep Field North (HDFN). It peaks 
around z ~ 0.8 in agreement with studies in other fields. 
The rest-frame 15 /im and bolometric infrared (8-1000 /im) 
luminosities of ISOCAM galaxies are computed using the 
correlations that we establish between the 6.75, 12, 15 /im 
and infrared (IR) luminosities of local galaxies. The result- 
ing IR luminosities were double-checked using radio (1.4 
GHz) flux densities from the ultra-deep VLA and WSRT 
surveys of the HDFN on a sample of 24 galaxies as well 
as on a sample of 109 local galaxies in common between 
ISOCAM and the NRAO VLA Sky Survey (NVSS). This 
comparison shows for the first time that MIR and radio 
luminosities correlate up to z ~ 1. This result validates the 
bolometric IR luminosities derived from MIR luminosities 
unless both the radio-far infrared (FIR) and the MIR-FIR 
correlations become invalid around z ~ 1. 

The fraction of IR light produced by active nuclei was 
computed from the cross-correlation with the deepest X- 
ray surveys from the Chandra and XMM-Newton obser- 
vatories in the HDFN and Lockman Hole respectively. We 
find that at most 20% of the 15 /im integrated galaxy 
light is due to active galactic nuclei (AGNs) unless a large 
population of AGNs was missed by Chandra and XMM- 
Newton. 

About 75 % of the ISOCAM galaxies are found to be- 
long to the class of luminous infrared galaxies (Lir > 10 11 
Lq). They exhibit star formation rates of the order of ~ 
100 Mq yr^ 1 . The comoving density of infrared light due 
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to these luminous IR galaxies was more than 40 times 
larger at z ~ 1 than today. 

The contribution of ISOCAM galaxies to the peak of 
the cosmic infrared background (CIRB) at 140 /im was 
computed from the MIR-FIR correlations for star form- 
ing galaxies and from the spectral energy distribution of 
the Seyfert 2, NGC 1068, for AGNs. We find that the 
galaxies unveiled by ISOCAM surveys are responsible for 
the bulk of the CIRB, i.e (16 ± 5) nW m~ 2 Hz" 1 as com- 
pared to the (25 ± 7) nW m~ 2 Hz -1 measured with the 
COBE satellite, with less than 10% due to AGNs. Since 
the CIRB contains most of the light radiated over the his- 
tory of star formation in the universe, this means that a 
large fraction of present-day stars must have formed dur- 
ing a dusty starburst event similar to those revealed by 
ISOCAM. 

Key words: Galaxies: evolution - Infrared: galaxies - 
Galaxies: starburst - Galaxies: Seyfert 



1. Introduction 

The extragalactic background light (EBL) is a measure- 
ment of the sum of the light produced by all extragalactic 
sources over cosmic time. When it is integrated over the 
full spectral range, the so-called cosmic background is a 
fossil record of the overall activity of all galaxies from 
their birth until now. It can be considered as the global 
energetic budget available for any model aiming at simu- 
lating the birth and fate of galaxies during the Hubble 
time. However the physical origin of this light will re- 
main unknown until we have pinpointed the individual 
sources responsible for it. The goal of the present paper 
is to demonstrate that an important new result has come 
from the combination of a series of deep extragalactic sur- 
veys performed in the mid-infrared (MIR) at 15 /im with 
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the ISOCAM camera (Cesarsky et al. 1996a) onboard the 
Infrared Space Observatory (ISO, Kessler et al. 1996): we 
suggest here that the galaxies detected in these surveys, 
which median redshift of z ~ 0.8 was measured from a 
sub-sample of ISOCAM galaxies, contribute dominantly 
to the cosmic infrared background (CIRB), i.e. the EBL 
integrated over all wavelengths within A= 5 to 1000 /im. 

The CIRB was recently detected and measured thanks 
to the cosmic background explorer (COBE) instruments 
FIRAS (Far Infrared Absolute Spectrometer) and DIRBE 
(Diffuse Infrared Background Experiment) (Puget et al. 
1996, Fixsen et al. 1998, Lagache et al. 1999, 2000, Hauser 
et al. 1998, Dwek et al. 1998, Finkbeiner et al. 2000) from 
100 /urn to 1 mm. It peaks around X max ~ 140 /jm and 
was found to represent at least half and maybe two thirds 
of the overall cosmic background (see Gispert, Lagache & 
Puget 2000). Hence the CIRB reflects the bulk of the star 
formation that took place over the history of the universe. 
By resolving it into individual galaxies, we would there- 
fore pinpoint the times and places where most stars seen 
in the local universe were formed. Two physical processes 
were considered for its origin: nucleosynthesis, i.e. stellar 
radiation in star forming galaxies, and accretion around 
a black hole, i.e. active galactic nuclei. In both cases, the 
light is not directly coming from its physical source but 
is reprocessed by dust, i.e. absorbed and re-radiated ther- 
mally by the "warm" dust. Both processes are probably 
related (see Genzel et al. 1998), but energetic considera- 
tions, based on the presence of massive black holes and on 
the amount of heavy elements in local galaxies, suggest 
that star formation should by far dominate in the CIRB 
over AGN activity (Madau & Pozzetti 2000, Franceschini 
et al. 2001). However, until the individual galaxies respon- 
sible for the CIRB are found and studied in detail, this 
result will remain theoretical. 

The spectral energy distribution (SED) in the IR of 
local galaxies peaks above ^60 /im and typically around 
80 ± 20 ^m (see Sanders & Mirabel 1996). As a result, 
the distant galaxies responsible for the peak of the CIRB 
detected by COBE around A max ~ 140 /im should be lo- 
cated below z rsj 1.3 and present a redshift distribution 
peaked around z ~ 0.8, if their SEDs do not strongly dif- 
fer from those of local galaxies. As we will see this is also 
the redshift range of the galaxies detected at 15 /tm with 
ISOCAM. 

The ISOCAM extragalactic surveys were performed 
with two filters, LW2 (5-8.5 /im) and LW3 (12-18 /im), 
centered at 6.75 and 15 /im respectively. The 6.75 /im sam- 
ple of sources is strongly contaminated by galactic stars, 
whereas stars are rather easily distinguished from galax- 
ies at 15 /im using optical-MIR colour-colour plots. As a 
consequence, we are only concerned here by the 15 /im 
galaxies. Moreover, the observed 6.75 /im light is no more 
produced by dust emission for galaxies more distant than 
z ~ 0.4 because of k-correction (redshifted stellar light 
dominates the 6.75 /im band above this redshift), whereas 



the observed 15 /im light is mostly due to dust emission 
for galaxies up to z ~ 2. 

About 1000 galaxies detected in the 15 /im surveys 
were used to produce number counts (i.e. surface density 
of galaxies as a function of flux density; see Elbaz et al. 
1999). The steep slope of the 15 /im counts below ~ 1 mJy 
indicates the presence of an excess of faint sources by one 
order of magnitude in comparison with predictions assum- 
ing no evolution of the 15 /im luminosity function with 
redshift. The presence of broad emission features in the 
MIR spectrum of galaxies alone cannot explain the shape 
of the number counts and a strong evolution of either the 
whole luminosity function (Xu 2000, Chary & Elbaz 2001) 
or preferentially of a sub-population of starburst galaxies 
evolving both in luminosity and density (Franceschini et 
al. 2001, Chary & Elbaz 2001, Xu et al. 2001) is required 
in order to fit the ISOCAM 15 /im counts. In the present 
paper, we suggest that these ISOCAM galaxies are in fact 
dusty starbursts responsible for the bulk of the CIRB. 

In Sect. |^, we compare the sensitivity of different ex- 
tragalactic surveys in several wavelength ranges to detect 
the galaxies responsible for the CIRB. It is suggested that 
MIR is presently the most efficient technique to detect 
dusty starbursts up to z ~ 1.3. 

In Sect. |, we calculate the 15 /im integrated galaxy 
light (IGL) due to ISOCAM galaxies. The 15 /im IGL is 
the sum of the 15 /im fluxes from individual galaxies, down 
to a given sensitivity limit, per unit area. It represents a 
lower limit to the 15 /im EBL, which remains unknown. 
Once the redshift distribution and SED of these galaxies 
is determined, it becomes possible to estimate their con- 
tribution to the CIRB. 

In Sect. |, we demonstrate that MIR luminosities at 
6.75, 12 and 15 /im are strongly correlated with the bolo- 
metric IR luminosity (from 8 to 1000 /im) for local galax- 
ies. The correlations presented in Chary & Elbaz (2001) 
are confirmed here with a larger sample of galaxies. 

Before spectroscopic redshifts are obtained for the full 
sample of ISOCAM galaxies used to produce these num- 
ber counts, the redshift distribution of these galaxies can 
be inferred from a few sub-samples: HDFN (Aussel et al. 
1999, 2001), CFRS-14 (Flores et al. 1999, 2002), CFRS- 
03 (Flores et al. 2002). The ultra-deep ISOCAM survey of 
the HDFN samples a flux density range where most of the 
evolution observed in the number counts takes place and 
where the bulk of the 15 /im IGL is produced. This field 
is complete in spectroscopic redshifts, so it is used to es- 
timate the bolometric IR luminosities and star formation 
rates of the ISOCAM galaxies in Sect. [|. 

This result relies on two assumptions: 

- that the main source for the MIR light in ISOCAM galax- 
ies is star formation and not accretion around a black hole. 

- that the correlations found in the local universe between 
the MIR and bolometric IR luminosity of galaxies remain 
valid up to z ~ 1. 
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The first assumption is discussed and justified in 2.2. Contribution of ISOPHOT galaxies 



Sect. 5.2, where soft and hard X-ray data from the Chan- 



dra and XMM-Newton X-ray observatories are combined 
with ISOCAM data on galaxies in the HDFN and Lock- 
man Hole regions respectively. 

The issue of the robustness of the MIR-FIR correla- 
tions in the distant universe is addressed in Sect. 5.5, 
where IR luminosities are also computed from radio (1.4 
GHz) flux densities for a sub-sample of 24 distant and 109 
local ISOCAM galaxies. 

In Sect. |5.6| , we compute the cosmic density of IR light 
due to luminous IR galaxies (£ir > 10 11 L®) at z ~ 1. 
In Sect. |, we evaluate the contribution of the ISOCAM 
galaxies to the CIRB, more precisely to its peak emission 
around X max ~ 140 fim. Finally, the nature of ISOCAM 
galaxies is discussed in the conclusions (Sect. 0). 

In the following, we will use the terms ULIG for galax- 
ies with an IR luminosity Lir = L[8 — 1000 /im] > 
10 12 L®, LIG, when 10 11 < (L m /L®) < 10 12 and lu- 
minous IR galaxies for both (Lir > 10 11 L®). Through- 



out this paper, we will assume H Q = 75 km 



Mpc~ 



a, 



,.= 0.3 and Q A = 0.7. 



2. How to resolve the CIRB into individual 
galaxies ? 

2.1. Contribution of IRAS galaxies 

One of the major results from the IRAS satellite was to 
demonstrate that the bolometric luminosity of galaxies 
more luminous than Lboi ~ 10 11 L® was underestimated 
by more than one order of magnitude before the IR lu- 
minosity was accounted for (Soifer et al. 1987). These 
luminous IR galaxies are not typical of the local galaxy 
population, since local galaxies radiate only about 30 % of 
their bolometric luminosity in the dust regime from 8 to 
1000 /im and only ~ 2 % of the local bolometric luminosity 
density is due to luminous IR galaxies. If we convert the 
IR luminosity into a star formation rate (SFR) using the 
formula of Kennicutt (1998), 

SFR(M Q /yr) = 1.71 x 10~ 10 L m [8 - 1000/xm](L Q ) (1) 

we find that luminous IR galaxies (Lt>oi ~ 
-^m > 10 11 L®) form stars at a rate larger than 
^20A/q yr _1 . Eq. ([!]) assumes continuous bursts last- 
ing 10-100 Myr, solar abundance and a Salpeter IMF 
(Kennicutt 1998). 

IRAS surveyed about 95 % of the sky down to a com- 
pleteness limit of ~ 0.5 mJy at 60 /im (~ 1.5 Jy at 100 /im) 
and the differential counts are well fitted by an Euclidean 
slope, dN/dS v ~ S- 2 - 5 (Soifer et al. 1987). The fluxes 
of all galaxies detected down to this sensitivity limit by 
IRAS add up to a 60 /mi IGL of - 0.15 nW m" 2 sr" 1 . 
This is less than 1 % of the value of the CIRB measured 



by COBE at A r , 



140 /im, IGL 140 = (25±7) nW m" 



Deeper FIR extragalactic surveys were performed with 
ISOPHOT (Lemke et al. 1996) onboard ISO at 170 /im. 
At this depth, the differential counts are no more fitted by 
an Euclidean slope, expected in the case of no evolution of 
galaxies with redshift, but instead a strong excess of faint 
sources was found. At the confusion limit of ISOPHOT of 
about 120 mJy, less than 10 % of the value of the CIRB 
at 170 /im is resolved into individual galaxies (Dole et al. 
2001). A preliminary follow-up of the ISOPHOT galaxies 
suggests that a large fraction of these galaxies are located 
at low redshift. 

2.3. Contribution of SCUBA galaxies 

In the sub-millimeter range, the SCUBA bolometer array 
on the James Clerk Maxwell Telescope (JCMT) has pro- 
duced deep images at 450 and 850 /im. Number counts at 
both wavelengths present a steep slope compatible with 
a strong evolution as compared to the local universe (see 
Smail et al. 2001 and references therein). 

At 450 /im, a depth of 10 mJy is reached and the com- 
bined fluxes of all SCUBA galaxies produce about 15 % of 
the CIRB measured by COBE-FIRAS at this wavelength 
(Smail et al. 2001). 

At 850 /^m, SCUBA is confusion limited at ~ 2 mJy 
(Hughes et al. 1998, Barger, Cowie & Sanders 1999, Eales 
et al. 2000, Smail et al. 2001), because of its large point 
spread function (PSF) of 15" full width half maximum 
(FWHM). About 20% of the value of the CIRB measured 
by COBE-FIRAS at 850 /im (the 850 /im EBL) is resolved 
into galaxies at this depth. However, using gravitational 
lensing this limit can be lowered to ~ 1 mJy, where ~ 60 % 
of the CIRB is resolved (Smail et al. 2001, Blain et al. 
1999a). 

The 850 /im EBL measured by COBE-FIRAS (~ 
0.5±0.2 nW m~ 2 sr -1 ) is 50 times lower than the peak 
value of the CIRB at X max ~ 140 /im. In order to com- 
pute the contribution of SCUBA galaxies to the peak of 
the CIRB it is therefore necessary to determine their red- 
shift distribution and SED. Until now, very few redshifts 
have been obtained due to the large PSF of SCUBA and 
optical faintness of these galaxies. Hence the fraction of 
the CIRB resolved into individual galaxies by SCUBA re- 
mains highly uncertain. 

2.4. Comparison with the MIR and radio 

The deepest ISOCAM surveys reach a completeness limit 
of ~ 100 /iJy at 15 /im (without lensing, but this limit 
goes down to ~ 30-50 /iJy including lensing, see Sect. ||). 
For a given galaxy SED and redshift, this flux den- 
sity can be converted into a bolometric IR luminosity 
(iiR [8-1000 /im]). Any galaxy more luminous than this 
Lir will be detected in the survey. In order to compare 
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A (/im) A (/im) 



Fig. 1. Spectral energy distribution of galaxies in the infrared, a) SEDs of M 82 (dashed line) and Arp 220 (plain 
line) combining photometric (filled and open dots) and spectroscopic data. Open dots are used when they overlap with 
spectroscopic data. Radio data for Arp 220 from Anantharamaiah et al. (2000). Sub-millimetric data for M 82 from 
SCUBA (Hughes et al. 1994). A = 45-190 /xm: ISO-LWS spectra from Fischer et al. (1999). Arp220's 18 /im silicate 
absorption feature from Smith, Aitken & Roche (1989). A = 5-18 /mi: ISOCAM Circular Variable Filter (CVF) spectra 
(Arp 220: Charmandaris et al. 1999, M 82: Forster Schreiber et al. 2000) normalized to the IRAS 12 /mi flux. For 
M 82, the CVF spectrum of the central region of the galaxy was complemented outward by adding to it a typical 
spectrum of a reflection nebula (NGC 7023, Cesarsky et al. 1996b), in order to match the 12 /im IRAS flux. A =2.5 
to 4/im: for M 82, SWS spectrum of Forster Schreiber et al. (2001) scaled using the NGC7023 spectrum of Moutou 
et al. (2000); for Arp 220, fit to the ISOPHOT data (Klaas et al. 1997) combined with the profile from Moutou et 
al. (2000). Optical to near-infrared: mix of the Coleman, Wu and Weedman (1980) templates to match the observed 
total magnitudes obtained from NED (Arp 220: 67% Sbc and 33% Ell; M 82: 12% Sbc and 88% Ell), b) 10 SEDs 
with log 10 (L m /L Q )= 8.5 to 13, with a step of 0.5, from the library of 105 SEDs of Chary & Elbaz (2001). 




Fig. 2. IR luminosity (left axis) and SFR (right axis) corresponding to the sensitivity limits of ISOCAM (15 /mi, plain 
line), ISOPHOT (170 /im, dotted line), SCUBA (850 /im, dot-dashed line) and VLA/WSRT (21 cm, dashed line) as 
a function of redshift. K-correction from: a) the template SED of M 82 (Fig. [l]a), normalized to the appropriate IR 
luminosity, b) the library of template SEDs from Chary & Elbaz (2001), plotted in the Fig. [l]b. 



the sensitivity to distant luminous IR galaxies of ISO- 
CAM, ISOPHOT and SCUBA, we have calculated L m 
corresponding to their sensitivity limits for a given red- 
shift between z=0 and 3. On one hand, we used the SED 



of M 82 (Fig. |T]a) for galaxies of all luminosities. On the 
other hand, we used a library of template SEDs from 
Chary & Elbaz (2001). In the following, we will call it 
the "multi-template" technique. This library of 105 SEDs 
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was constructed under the constraint that it reproduces 
the correlations observed between the 6.75, 12, 15, 60, 
100 and 850 fim luminosities of local galaxies. The SEDs 
were interpolated from a sample of observed SEDs, includ- 
ing MIR spectra obtained with the ISOCAM CVF. They 
cover the luminosity range logio (Lir/Lq) = 8.5 to 13.5 
(see Fig. 0b). 

We did not use Arp 220's SED in Fig. 0a, because this 
galaxy presents a large FIR over MIR luminosity ratio, 
which is not representative of galaxies within the same 
luminosity range (see Figs.||c,d). 

The radio continuum (A= 21 cm, v= 1.4 GHz) is also 
a tracer of L\r because of the tight correlation between 
both luminosities in the local universe (see Yun, Reddy 
& Condon 2001 and references therein). If this correla- 
tion remains valid in the distant universe then we can 
translate the sensitivity limit of the deepest radio surveys 
into a minimum Lir accessible at a given redshift assum- 
ing a radio SED. The correlation between the radio and 
FIR luminosities is usually described by the "g" parameter 
(Condon et al. 1991): 

, ( £fir(W) 1 \ ro s 

q - 1 ° 9W 1.3.75 x 10-(Hz) X L.^fWHr 1 )] (2) 

where Li.4GHz is the monochromatic luminosity at 1.4 
GHz and Lfir is the FIR luminosity between 40 and 
120 /mi, as defined by Helou et al. (1988), 

L F1R = 1.26 x 10" 14 (2.58 S 60 + S WQ ) x 47rd(m) 2 (3) 

where Seo and «Sioo are the IRAS flux densities at 60 and 
100 /im in Jy. The relationship between Lfir [40- 120 /mi] 
and Lir[8-1000 /im] was computed with the sample of 300 
galaxies from the Bright Galaxy Sample (BGS, Soifer et 
al. 1987) detected in the four IRAS bands (12, 25, 60 and 
100 /mi): 

L m = (1.91 ± 0.17) x L FIR (4) 

Lir is defined as (Sanders & Mirabel 1996): 

Lir = [13.48 5i2 + 5.16 S 25 + 2.58 S 60 + S wo ] (t ~ 
x 1.8 x 10~ 14 x 47rd(m) 2 ( ' 

where S12 and S25 are the IRAS flux densities at 12 and 
25 /mi in Jy. Yun, Reddy & Condon (2001) measured a 
value of q= 2.34±0.01 from a flux limited (S 60 > 2 Jy) 
sample of 1809 IRAS galaxies. 

The deepest radio surveys presently available were per- 
formed in the HDFN by Richards (2000) and Garrett et 
al. (2000) with the VLA and WSRT respectively. Both 
surveys reach the same sensitivity limit of 40/iJy at 1.4 
GHz (5-er). In order to convert this flux density into a 
-£a.4GHz for a given redshift, we used a power index of 
a— 0.8±0.15 (S„ (x ^~ Q ), as suggested by Yun, Reddy & 
Condon (2001) for starburst galaxies. 

Fig. H presents the sensitivity of the deepest MIR, FIR, 
sub-millimeter and radio surveys in the form of Lir (or 



SFR) as a function of redshift. The sensitivity limits used 
for ISOCAM, ISOPHOT and SCUBA are taken from the 
deepest existing surveys at these wavelengths, as previ- 
ously described. Fig. ||a (with M 82) and Fig. ("multi- 
template" technique) both clearly show that the faintest 
IR galaxies are best detected at 15 /mi up to z ~ 1.3. The 
right axis of the plots shows the corresponding minimum 
SFR that a galaxy must harbor to be detected at a given 
redshift in the surveys. At z ~ 1, ISOCAM (plain line) 
detects all LIGs while SCUBA (dot-dashed line) detects 
only ULIGs. The difference in sensitivity between ISO- 
CAM and the deepest radio surveys is about a factor 2. 
ISOPHOT (dotted line) is only sensitive to ULIGs above 
z ~ 0.5. In Fig. ||b, where we used the library of template 
SEDs from Chary & Elbaz (2001), the plain line corre- 
sponding to ISOCAM rises faster above z ~ 1 than in 
Fig. |a, where we used M 82's SED. This different behav- 
ior comes from the fact that at these redshifts ISOCAM 
measures fluxes at about 7 /im in the rest-frame of the 
galaxy and we will show that Lir/ 'L7 Mm increases with 
Lir (Sect. |). 

Above z~2, the sub-millimeter becomes the most ef- 
ficient technique, although only galaxies more luminous 
than a few 10 12 Lq can be detected above a sensitivity 
limit of 2 mJy at 850 /zm. Hence the unlensed ISOCAM 
and SCUBA surveys are not sampling the same redshift 
and luminosity ranges. This statement is confirmed obser- 
vationally in the HDFN itself where none of the ISOCAM 
sources are detected by SCUBA (Hughes et al. 1998) and 
over a larger scale in the Canada France Redshift Sur- 
vey 14 (CFRS-14, Eales et al. 2000). In this latter field of 
~50' 2 , only the two brightest 15 /im ISOCAM sources are 
detected at 850 /im, among a sample of 50 ISOCAM (Flo- 
res et al. 1999) and 19 SCUBA sources. This confirms that 
both instruments detect different sets of objects. More 
common objects between ISOCAM and SCUBA are ex- 
pected if they are gravitationally lensed, since the sensi- 
tivity limits of both instruments are then decreased by a 
factor of about two. Indeed in the clusters A 370 and A 
2390, three of the four lensed SCUBA galaxies (Smail et 
al. 2001) are also detected with ISOCAM (Altieri et al. 
1999, Metcalfe 2000). 

3. The 15 fim extragalactic background light 

Above the Earth's atmosphere, the 15 /mi light is strongly 
dominated by the zodiacal emission from interplanetary 
dust and it has not yet been possible to make a direct 
measurement of the 15 /im background, or EBL. Individ- 
ual galaxies contribute to this background and a lower 
limit to the 15 /mi EBL can be obtained by adding up 
the fluxes of all ISOCAM galaxies detected per unit area 
down to a given flux limit. The resulting value is called 
the 15 /xm integrated galaxy light (IGL). 

The statistical reliability with which the surface den- 
sity of MIR sources is computed depends on their flux: 
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Fig. 3. a) Differential contribution to the 15 /im Integrated Galaxy Light as a function of flux density and AB magni- 
tude. The plain line is a fit to the data: Abell 2390 (Altieri et al. 1999), the ISOCAM Guaranteed Time Extragalactic 
Surveys (IGTES, Elbaz et al. 1999), the European Large Area Infrared Survey (ELAIS, Serjeant et al. 2000) and 
the IRAS all sky survey (Rush, Malkan & Spinoglio 1993). b) Contribution of ISOCAM galaxies to the 15 /im ex- 
tragalactic background light (EBL), i.e. 15/im Integrated Galaxy Light (IGL), as a function of sensitivity or AB 
magnitude (AB = —2.5log(S nl jy) + 16.4). The plain line is the integral of the fit to dIGL/dS (Figj|a). The dashed 
lines correspond to 1-er error bars obtained by fitting the 1-cr upper and lower limits of dIGL/dS. 



large and shallow fields are required for bright sources 
which are less numerous, whereas faint sources require 
smaller and deeper fields. The optimization is obtained by 
observing similar volumes of the sky but in pencil beams 
for faint sources and in wide surveys for bright sources, 
instead of one large and deep survey which would be un- 
necessarily time consuming. Together with the need to 
account for cosmic variance, this explains why we have 
combined a set of 10 surveys located in 7 regions of the 
northern and southern hemispheres, in order to estimate 
the 15 /im IGL. These surveys cover five orders of magni- 
tude in flux density from ~ 50 /iJy to ~ 4 Jy (including 
IRAS counts at 12 /im). The number counts presented in 
Elbaz et al. (1999) are converted here into a differential 
contribution to the 15 /zm IGL as a function of flux den- 
sity, estimated from the following formula: 



dIGL 
dS 



dN 
~d~S 



(G) 



where dN(si *) is the surface density of sources with a 
flux density 5^ [15 /im]= S15 (mJy) over a flux density bin 
dS(mJy) (1 mJy= 10~ 20 nW m~ 2 Hz" 1 ) and is 15 (Hz) is 
the frequency of the 15 /im photons. 

Below S15 ~ 3 mJy, about 600 galaxies were used to 
produce the points with errors bars in Fig. ||a. Fig. ||b 
shows the 15 /im IGL as a function of depth. It corresponds 
to the integral of Fig. ||a, where the data below 3 mJy were 
fitted with a polynomial of degree 3 and the 1-cr error bars 
on dIGL/dS were obtained from the polynomial fit to 
the upper and lower limits of the data points. The 15 /im 
IGL does not converge above a sensitivity limit of S15 ~ 



50/iJy, but the flattening of the curve below S15 ~ 0.4 
mJy suggests that most of the 15 /jm EBL should arise 
from the galaxies already unveiled by ISOCAM. 

Above the completeness limit of S15 ~ 50 /iJy, we com- 
puted an IGL of: 



IGL 15 (S 15 > 50 /zJy) = 2.4 ± 0.5 nW m 



- 2 sr- 1 



(7) 



where the error bar corresponds to the 68 % confidence 
level (i.e. 1-cr). The error bar combines the uncertainty on 
the measurements for each individual survey with cosmic 
variance. Each flux bin is covered by two to three indepen- 
dent surveys. ISOCAM number counts are only complete 
above S15 ~ 100 /iJy, but Aussel et al. (1999) computed 
number counts down to S15 ~ 50 /iJy after correcting for 
incompleteness and a consistent result is obtained by Al- 
tieri et al. (1999) after correcting for lensing magnification 
a sample of galaxies in the field of the galaxy cluster Abell 
2390. 

Only 9 sources are detected below S15 ~ 50 /iJy by Al- 
tieri et al. (1999), where the flux determination and correc- 
tion for completeness are less robust and are not confirmed 
by another survey. A tentative value of IGL15 ~ 3.3 ±1.3 
nW m~ 2 sr -1 was quoted in Altieri et al. (1999), above 
Si5= 30/iJy. This value may be slightly overestimated if 
we consider that the models of Franceschini et al. (2001) 
and Chary & Elbaz (2001), which reproduce the number 
counts from ISOCAM at 15 /mi, from ISOPHOT at 90 and 
170 /im and from SCUBA at 850 /im, as well as the shape 
of the CIRB from 100 to 1000 /im, consistently predict a 
15 /im EBL of: 



EBL modeU (15fim) ~ 3.3 nW m~ 2 sr^ 1 



(8) 
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Field 


F(mJy) 


% 


Ref. 


A2390,HDFN 


[ 0.05, 0.1 ] 


23 


(1,2) 


IGTES 


[ 0.1, 0.5 ] 


48 


(2,3) 


IGTES 


[ 0.5, 3 ] 


23 


(3) 


ELAIS 


[ 3, 32 ] 


4.5 


(4) 


Interpolation 


[ 32,300 ] 


1 


(5) 


IRAS 


[ 300,4000 ] 


0.5 


(6) 



Table 1. Relative contributions to the 15 /im IGL as a 
function of flux density. Ref.: (1) Altieri et al. (1999), Met- 
calfe (2000), (2) Aussel et al. (1999, 2001), (3) Elbaz et al. 
(1999), (4) Serjeant et al. (2000), renormalized by a factor 
0.5 as in Genzel & Cesarsky (2000), following suggestion 
from S.Serjeant for photometric calibration, (5) Interpola- 
tion in a flux density range where we do not have enough 
statistics. We used the same slope as in the IRAS counts, 
(6) IRAS counts from Rush, Malkan & Spinoglio (1993) 
at 12 /im. 



If this prediction from the models is correct then about 
73 ± 15 % of the 15 /im EBL is resolved into individual 
galaxies by the ISOCAM surveys. 

This result is consistent with the upper limit on the 
15 fim EBL estimated by Stanev & Franceschini (1998) 
of: 



EBL ma3 : (15/zm) 



5 nW m 2 sr 1 



(9) 



This upper limit was computed from the 1997 7-ray out- 
burst of the blazar Mkn 501 (z = 0.034) as a result of 
the opacity of MIR photons to 7-ray photons, which an- 
nihilate with them through e + e _ pair production. It was 
recently confirmed by Renault et al. (2001), who found an 
upper limit of 4.7 nW m -2 sr^ 1 from 5 to 15 /im. 

Finally, the 15 /an background (EBL) must be con- 
tained between the following limits: 



2.4 ±0.5 < 



EBL(15nm) 



(nW m- 2 sr- 1 ) " 5 ^ 

In Table |, we give the relative contribution to the 
15 /im IGL from sources within a given flux density range. 
Nearly all the IGL is produced by sources fainter than 3 
mJy (94%) and about 70% by sources fainter than 0.5 
mJy. This means that the nature and redshift distribu- 
tion of the galaxies producing the bulk of the 15 /im IGL 
can be determined by studying these faint galaxies only. 
We will therefore assume in the following that the spectro- 
scopically complete sample of 15 /jm galaxies detected in 
the HDFN below 0.5 mJy (Aussel et al. 1999, 2001) is rep- 
resentative of the population of dusty galaxies producing 
the bulk of the 15 /im IGL. 



4. MIR, as a tracer of the bolometric IR 
luminosity 
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Fig. 4. Spectral energy distributions in the infrared (5- 
200 /im) of 3 starbursts (M 82, Arp 244, Arp 220) and 
2 Seyfert 2's (NGC 1068, NGC 6240). MIR spectra from 
ISOCAM CVF: NGC 1068 (Laurent et al. 2000, Le Floc'h 
et al. 2001), Arp 244 (Vigroux et al. 1996, Mirabel et al. 
1998), NGC 6240 (Charmandaris et al. 1999). FIR: fit of 
the IRAS flux densities at 12, 25, 60 and 100 /im using the 
template SEDs from Chary & Elbaz (2001) for Arp 244 
and NGC 6240. For NGC 1068, the IRAS flux densities 
were fitted by the sum of two modifified black bodies of 
temperatures 30 and 170 K as in Telesco et al. (1984). 
Open dots with error bars: IRAS data (from NED) for 
Arp244, NGC 6240 and NGC 1068. 



4-1- SED of galaxies in the IR 

Above a bolometric luminosity of a few 10 10 Lq, the SED 
of a galaxy is dominated by dust radiation with respect 
to direct stellar emission. Stellar light is absorbed by dust 
and re-radiated in the IR regime above A ~ 5 /im. The full 
SEDs, from the UV to the radio, of the starburst galaxy, 
M 82, and of the ULIG, Arp 220, are plotted in the Fig. fj]a. 
A zoom on the 5-2000 /im wavelength range of the SEDs 
is shown in the FigJJ, together with three other galaxies: a 
starburst, Arp 244, and two Seyfert 2 galaxies, NGC 1068 
and NGC 6240. Arp 244, also known as the "Antennae" 
galaxy, is a system of two interacting spiral galaxies (Vi- 
groux et al. 1996, Mirabel et al. 1998). The optical, MIR 
and FIR luminosities of these five galaxies are given in 
the Table ||. Note the increase of the MIR luminosity to- 
gether with i F iR, or L m , from M 82 to Arp 244 to NGC 
6240. This correlation will be statistically established in 



Sect. 4.2.2 



The FIR over MIR luminosity ratio of Arp 220 is larger 
than for galaxies of equivalent luminosity (see Fig. [5] in 
Sect. p~2~2; ). The flat FIR over MIR ratio for NGC 1068 
is typical of the hot dust continuum from AGNs. NGC 
1068 is the closest Seyfert 2 (L m = 1.7 x 10 11 Lq). On 
one hand, its MIR spectrum is completely dominated by a 
continuum due to the hot dust heated by the central AGN, 
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M 82 


Arp244 


Arp220 


N1068 


N6240 


d(Mpc) 


3.25 


22.0 


78.4 


14.4 


99.6 


log(L B ) 


10.2 


10.2 


10.6 


10.5 


10.5 


log(Lm) 


10.6 


10.8 


12.2 


11.2 


11.8 


log{L F iR) 


10.4 


10.7 


12.1 


10.8 


11.7 


log(L LW s) 


9.8 


10.0 


10.6 


10.8 


10.7 


log(L L w2) 


9.9 


10.1 


10.3 


10.6 


10.5 


LmJ Lb 


2 


3 


39 


5 


18 


Lm/ Llwz 


15 


7 


41 


5 


10 


LmJ Llw2 


4 


4 


87 


4 


22 



Table 2. Luminosities of 3 starbursts (M 82, Arp 244, 
Arp 220) and 2 Seyfert 2's (NGC 1068, NGC 6240). L B is 
the B-band monochormatic luminosity {yL v at 0.44 /im) 
from RC3 (de Vaucouleurs et al. 1991). Lir and Lfir 
are the IR and FIR luminosities integrated over 8-1000 
and 40-500 /im respectively from the spectra shown in the 
Fig.^. Llw3 and Lj_,w2 are the monochromatic luminosi- 
ties {vL v ) at the central wavelengths 15 and 6.75 /xm of 
the ISOCAM LW3 and LW2 filters respectively. At z~ 
1, the LW3 filter (12-18 /im) matches the rest- frame LW2 
bandpass (5-8.5 /im). 



which produces about 75 % of the MIR luminosity (Le 
Floc'h et al. 2001). On the other hand, the bulk of its FIR 
luminosity originates from the diffuse region surrounding 
the nucleus as shown by the SCUBA image at 450 /im, 
which is associated with star formation (Le Floc'h et al. 
2001). As a result, the FIR over MIR luminosity ratio is 
much lower for NGC 1068 than for starbursts (see Fig.|J). 
This galaxy is a perfect illustration of the close connection 
between star formation and AGN activity. 

The MIR (5-40 /im) and FIR (40-1000 /im) emission of 
a galaxy combine three major components (see Fig^): 

- broad emission features and their associated under- 
lying continuum, which dominate in the A ~ 5-10 /im do- 
main, and up to ~ 20 /im for galaxies with moderate star 
formation. These bands located at 6.2, 7.7, 8.6, 11.3 and 
12.7 /im, are usually quoted as PAHs (polycyclic aromatic 
hydrocarbons, Leger & Puget 1984, Puget & Leger 1989, 
Allamandola et al. 1989), although their exact nature re- 
mains uncertain, e.g. Jones & d'Hendecourt (2000). 

- the "warm" dust continuum produced by very small 
dust grains (VSGs, smaller than 0.01 /im, Andriesse 1978, 
Desert, Boulanger & Puget 1990) heated at temperatures 
larger than 200 K without reaching thermal equilibrium. 
This component dominates the SED between 10-20 to 
40 /im for luminous galaxies. 

- the "cold" dust continuum produced by big dust 
grains ("big grains", of size > 0.01 /im) in the thermal 
equilibrium at cool temperatures (below 60-70 K typi- 
cally) . This component is responsible for the bulk of the 
FIR light where the SED peaks. 

4.2. The MIR- FIR correlation 



4.2.1. The database 

We have compiled a database containing MIR and FIR 
fluxes for a sample of 157 galaxies detected with ISOCAM 
and IRAS: 50 galaxies from Boselli et al. (1998), 44 galax- 
ies from Roussel et al. (2001), 19 galaxies from Laurent 
et al. (1999, 2000), 41 galaxies from Dale et al. (2001), 
plus M 82, Arp 244 and NGC 1068 (from Figs. |). Fig. § 
shows that the 6.75, 12, 15 and bolometric IR (8-1000 /im) 
luminosities (£ir) of these galaxies are tightly correlated 
with each other. Lm is calculated following Eq. (||). At 
15 /im (LW3 filter), we completed our sample with the 
NEPR sample of IRAS galaxies observed with ISOCAM- 
LW3 (Aussel et al. 2000). We converted the 60 /im lumi- 
nosities of the NEPR galaxies into Lir by fitting the IR 
versus 60 /im correlation obtained with the BGS (Soifer et 
al. 1987). In order to complete our sample in the high lu- 
minosity range at 6.75 /im (LW2 filter), we used a sample 
of 40 ULIGs observed with the spectrometer of ISOPHOT, 
PHOT-S (from Rigopoulou et al. 1999: the conversion into 
LW2 flux densities was computed as in Chary & Elbaz 
2001). We also used four ULIGs observed with the ISO- 
CAM CVF from Tran et al. (2001). We conservatively 
selected these star forming galaxies with respect to AGNs 
from their line over continuum ratioas suggested in Tran 
et al. (2001). The ISOCAM-CVF spectra were converted 
into LW2 flux densities using the transmission curve of the 
LW2 filter. 

4.2.2. Correlations 

The 6.75, 12 and 15 /im luminosities of local galaxies are 
tightly correlated with each other (see Fig. ||a,b) following 
the formulae: 

iAL„[15/zm] = 0.43 x vL v [6.75 fj,m] 

for ^£^[6.75 /im] < 5 x 10 9 Lq . . 
= 3.8 x 10- 7 x (i/L„[6.75 /im]) 162 
for vL„[6.75 /im] > 5 x 1O 9 L 

vL v [15nm] = 0.042 x (^[12 /mi]) 1 - 12 (12) 

These correlations can be used to estimate the rest-frame 
15 /im luminosity for galaxies located around z ~ 0.3 and 
1. At these redshifts, 15 /im observed corresponds to about 
12 and 7 /im in the rest-frame. 

The MIR luminosity of galaxies is a good tracer of Lir 
as seen in Figs. |^c,d, where Lir is correlated with vL v at 
15 and 6.75 /im following the formulae: 

Lir = 11.1+13 x (^[15 /im]) - 998 (13) 

L m = 4.78 +2 ;gx (^,[6.75 /im]) - 998 

for i/L„[6.75/Ltm] < 5 x 10 9 Lq , 
= 4.37± 2 -gx HT 6 x (^[6.75 /im]) 162 1 ' 
for uL v [6.75 /im] > 5 x 10 9 L Q 

The formulae used in these four correlations are con- 
sistent with Chary & Elbaz (2001), although slightly more 
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vL„[15 /xm] (L ) 1/LJ6.75 fan] (L s ) 



Fig. 5. IR luminosity correlations for local galaxies (see text for the origin of the sample). The five galaxies from Figjj 
and Table || are marked with open squares, a) ISOCAM-LW3 (15 /im) versus ISOCAM-LW2 (6.75 /im) luminosities 
iyL v ) (56 galaxies), b) ISOCAM-LW3 (15 /im) versus IRAS-12 /mi luminosities (45 galaxies), c) Lir[8-1000 /mi] versus 
ISOCAM-LW3 (15 /mi) luminosity (120 galaxies), d) LiR,[8-1000/mi] versus LW2-6.75/mi luminosities (91 galaxies). 
Filled dots: galaxies from the ISOCAM guaranteed time (47 galaxies including the open squares). Open dots: 40 galaxies 
from Rigopoulou et al. (1999). Empty triangles: 4 galaxies from Tran et al. (2001). Galaxies below Lir ~ 10 10 Lq 
present a flatter slope and have Lj-r/Lb < !• 



data were used here. The change of slope in Fig. |5|d 
and Eq. ( |i"4| ) from a slope of 1 to a slope larger than 
1 statistically establishes what was suggested by Fig. |]: 
Lir./^L„[6.75 /mi] increases with Lir, for galaxies with 
L m > 2xl0 10 Lq. 

The position of the five templates from Fig. [| are indi- 
cated with open squares on Fig. Two galaxies do not fit 
the correlations: Arp 220 and NGC 1068. Arp 220 presents 
a very large FIR over MIR ratio, which was interpreted as 
being due to extinction in the MIR by Haas et al. (2001), 
while the IR SED of the Seyfert 2 NGC 1068 is fiat, which 
is typical of AGNs. In this sense NGC 6240 is quite atyp- 
ical, since it falls perfectly on the correlation followed by 
star forming galaxies. The presence of an active nucleus 
in NGC 6240 was revealed by its strong X-ray luminos- 
ity and Ikebe et al. (2000) suggested that a large fraction 



of the IR luminosity of this galaxy could originate from 
its AGN. However, the presence of strong MIR features 
(Genzel et al. 1998), the absence of a strong continuum 
emission below A ~ 6/zm (Laurent et al. 2000), together 
with a FIR over MIR ratio typical of star forming galax- 
ies, suggest that the bulk of its IR luminosity may still be 
due to star formation. 

5. IR luminosities and star formation rates of 
ISOCAM galaxies 

In this Section, we compute the IR luminosity and star for- 
mation rate of the galaxies with S15 > 0.1 mJy located in 
the Hubble Deep Field North plus its flanking fields, here- 
after called ISOCAM-HDFN galaxies (Aussel et al. 1999, 
2001). This field presents the advantage of being spectro- 
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Fig. 6. a) Histogram of the flux densities of the 41 ISOCAM-HDFN galaxies (£15 > 0.1 mJy, Aussel et al. 1999, 2001) 
with poissonian error bars. Only one galaxy has no spectroscopic redshift (dashed line), b) Redshift distribution of 
the 40 galaxies with spectroscopic redshift (98% complete in redshift, poissonian error bars). 



scopically complete and deep enough to cover a flux den- 
sity range (S15 < 0.5 mJy) where more than 70% of the 
15 fim IGL is produced. Moreover, the redshift distribu- 
tion of the ISOCAM-HDFN galaxies (z ~ 0.8) is similar 
to the one measured in the CFRS fields CFRS-14 (Flo- 
res et al. 1999, 2002) and CFRS-03 (Flores et al. 2002). 
Together, these fields of 100' 2 each, possess 81 galaxies 
with Sis > 0.35 mJy. The redshift distribution of 70% of 
the CFRS galaxies, for which a spectroscopic redshift was 
measured, is very similar to the one presented below for 
the ISOCAM-HDFN galaxies and peaks around z ~ 0.7. 

5.1. Redshift distribution of the ISOCAM-HDFN galaxies 

The list of ISOCAM-HDFN sources used in the following 
comes from a revised version of the list published in Aussel 
et al. (1999). In this version (Aussel et al. 2001), 4 ISO- 
CAM sources were deblended using PSF deconvolution 
based on the position of the optical counterparts and the 
flux density of 10 other galaxies were re-evaluated. With 
a PSF FWHM of 4.6" at Ibfxm (Okumura 1998), opti- 
cal counterparts are easily identified for ISOCAM galaxies 
and all of them are found within 3" of the PSF centroid. 
This can also be seen is the crowded field of the nearby 
galaxy cluster Abell 1689 (Fadda et al. 2000). More than 
90 % of the ISOCAM galaxies present an optical counter- 
part brighter than I(Kron-Cousins)=22.5, which is natu- 
ral since ISOCAM is mostly sensitive to galaxies closer 
than z ~ 1.3 (see Fig. ||). It results from a combination of 
sensitivity limit and k-correction: dust emission becomes 
negligible below ~ 6 /im (except for bright AGNs which 
can be detected further away). 

A total of 41 ISOCAM-HDFN galaxies have 15 ^m flux 
densities between S15 =0.1 and 0.5 mJy (see Fig.^a) 
over a field of 25.8' 2 (HDFN+FF). Monte-Carlo simu- 



lations give a completeness of ~ 90 % in this flux den- 
sity range (Aussel et al. 1999, 2001). All but one galaxy 
(HDF_PS3_6a) possess a spectroscopic redshift (98 % com- 
plete). Their mean and median redshift is z ~ 0.8. Except 
for one galaxy located at z=2.01 (HDF_PM3_5, classified 
as an AGN, see Sect. [T^) and four galaxies closer than 
z=0.15, all ISOCAM galaxies have a redshift between z ~ 
0.4 and 1.3 (see FigJ|b), as expected from the Fig. |[ 

5.2. Fraction of AGN galaxies 

Cohen et al. (2000) classified the optical spectra of 2/3 of 
the ISOCAM galaxies and found only two galaxies pre- 
senting broad emission lines typical of AGNs: 
HDF.PM3_6 (C36367.1346, z = 0.960), 
HDF.PM3_20 (C36463_1404, z = 0.962). 
However, dusty AGNs may not be easily identified in 
the optical. 

The best technique at present for the identification of 
dusty AGNs is offered by hard X-ray observations, which 
are much less affected by extinction. The deepest soft plus 
hard X-ray observation currently available was performed 
with the Chandra X-ray Observatory in the HDFN by 
Hornschemeier et al. (2001) and Brandt et al. (2001,2002). 
The Chandra and ISOCAM catalogs were cross-correlated 
by Fadda et al. (2001). A total of 16 ISOCAM-HDFN 
(S15 > 0.1 mJy) are detected in the 1 Msec Chandra im- 
age at the depth of 5[0.5-2 keV]~ 3xl0~ 17 erg cm* 2 s" 1 
and S*[2-8 keV]~ 2xl0~ 16 erg cm" 2 s -1 . They represent 
39% of the ISOCAM-HDFN galaxies. However, the Chan- 
dra image is so deep that it is able to detect starburst 
galaxies up to large redshifts and among the 16 ISOCAM 
galaxies in common with Chandra sources, only 5 are clas- 
sified as being AGNs or AGN dominated by Fadda et al. 
(2001). These 5 sources include the two galaxies previ- 
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Fig. 7. Rest-frame MIR luminosity of ISOCAM-HDFN galaxies, a) luminosity of ISOCAM galaxies at A rf = 
15/(l+z) /im as a function of A r f. The clotted lines represent the redshifted LW3 (12-18 /mi) band, b) rest-frame 
15/im luminosity {vL v ) versus redshift estimated with three different techniques: M 82 SED (crosses), "closest filter" 
(filled dots) and "multi-template" (open dots). 



ously identified as AGNs from their optical spectra and 
the following 3 galaxies: 

HDF.PM3.5 (J123635.6+621424, z = 2.010), 

HDF.PS3.10 (J123642. 1+621545, z = 0.857), 

HDF.PM3.21 (J123646.4+621529, z = 0.851). 

Galaxies detected in the ultra-hard X-ray band (4-8 
keV) or with Lx > 10 43 erg s _1 were classified as AGNs. 
Galaxies with Lx < 10 40 erg s _1 or which could be fitted 
with the SED of M 82 or Arp 220 were classified in the 
starburst category. 

Finally only (12±5) % (5/41 galaxies, poissonian er- 
ror bar) of the ISOCAM-HDFN galaxies are found to be 
powered by an active nucleus. These AGNs contribute to 
(18±7) % of the sum of the flux densities (the IGL) of 
all 41 ISOCAM-HDFN galaxies. Since the active nuclei in 
these galaxies may not be responsible for the totality of 
the IR emission, this result suggests that at most (18±7) % 
of the 15/im IGL due to galaxies with Sis= 0.1-0.5 mJy 
is produced by AGN activity. 

The Si5= 0.5-3 mJy flux density range is covered by 
the ISOCAM image of the Lockman Hole, where Fadda 
et al. (2001) cross-correlated ISOCAM with the XMM- 
Newton X-ray observatory catalog of sources brighter than 
S[2-10 keV]~ 1.4xl0" 15 erg cm" 2 s" 1 and S[5-10 keV]~ 
2.4xl0 -15 erg cm -2 s _1 (Hasinger et al. 2001). Among 
103 ISOCAM sources with a flux density ranging from 0.5 
to 3 mJy, 13 sources were found to be AGN dominated 
by Fadda et al. (2001), i.e. (13±4)%. Such galaxies could 
be fitted by the SEDs of either type-I or type-II AGNs 
by Franceschini et al. (2002). These AGNs contribute to 
(15±5) % of the 15 /im IGL due to galaxies in the 0.5-3 
mJy flux range. 

The combination of the HDFN and Lockman Hole 
cover a flux density range of 0.1-3 mJy where 70% of 



the 15 /im IGL is measured. In this range, these two fields 
suggest that AGNs contribute to at most (17±6)%. It 
must be noted that the XMM-Newton and Chandra ex- 
periments were limited to energies below 10 keV. Some 
harder X-ray AGNs may not have been detected by these 
instruments yet and may be common with our ISOCAM 
source list. The presence of such objects will remain will 
highly uncertain until the launch of the next generation 
X-ray experiments (XEUS, Constellation X). We will as- 
sume in the following that the AGN fraction computed 
in the HDFN (which is larger than in the Lockman Hole) 
applyies for the 15 /im IGL as a whole. 

5.3. Rest-frame 15 [im luminosity of ISOCAM galaxies 

In this section, we compare three different techniques 
to estimate the rest-frame 15 /im luminosity of ISOCAM 
galaxies. The first and simplest technique consists in us- 
ing a single template SED to compute the k-correction. 
We chose the SED of M 82 (Fig. |l]a) which falls on the 
correlations plotted in the Fig. || (crosses in Fig. 0b). As 
in Sect, g, we decided not to use the SED of Arp 220 
(Fig. |l]a) because it is well outside the correlations found 
for galaxies of similar luminosities with its large FIR over 
MIR luminosity ratio. 

In a second technique, that we will call "closest fil- 
ter", we use directly the correlations described in Sect. |^ 
(Eq. ( pl| ) and Eq. (p^)) between the rest-frame 15 /im and 
6.75 or 12 /im luminosities. We first calculate the rest- 
frame wavelength A r f = 15/(1+2; ) /im and find which of the 
three filters ISOCAM-LW3, ISOCAM-LW2 and IRAS-12 
has its central wavelength the closest to A r f (see F ig, ^a). 
Then we apply the formula given in the Sect. |4.2.2| for this 
"closest filter" to compute the rcst-framc 15 /im luminos- 



12 



D. Elbaz et al 



Table 3. Properties of the ISOCAM galaxies detected above 0.1 mJy in the HDFN+FF (40 galaxies with spectroscopic 
redshift within a total of 41 galaxies). Col.(l): the galaxies have been divided into 5 categories (see text). Col. (2): 
median redshift. Col. (3): logarithm of the median 15 /im luminosity (£15 = vL v [15 /im]). Col. (4): logarithm of the 
median IR luminosity (8-1000 /im), calculated from the multi-template technique. Col. (5): median star formation rate, 
calculated from Lir using Eq. (|l|). The l-cr error bars in Cols. (3), (4), (5) are the median of the l-cr error bars for 
each individual galaxy. They characterize the uncertainty on L15, Ltr and (SFR) for individual galaxies, not their 
dispersion around the median value. Col. (6): number of galaxies in each category. Col. (7): relative contribution of each 
galaxy type to the 15 /jm IGL (/GL15) produced by the sum of all 40 galaxies. The mean and median values given in 
the last two lines concern only the 35 star forming galaxies. 



Type 


2 


log{Li 5 ) 


log(Lm) 


(SFR) 


Nb 


% IGL15 






(Lq) 


(Lq) 


(M & /yr) 






ULIGs 
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Normal 
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8.3± 0.1 


9.4± 0.2 
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11 


AGNs 


0.9 


11. 1± 0.7 


11.5± 0.7 




5 


19 


Median 


0.8 


10.4± 0.1 


11.4± 0.2 


41±19 


35 


100 


Mean 


0.8 


10.8± 0.1 


11.8± 0.2 


109±51 


35 


100 



ity (filled dots in Fig. 0b). Since most ISOCAM-HDFN 
galaxies are located around the median redshift z ~ 0.8, 
the correlation between LW3 and LW2 (Fig. |d, Eq. (fl]) 
is most often used. 

The third techni que is the multi-template technique 
already used in Sect. 2^4, where we use the library of 105 
template SEDs from Chary & Elbaz (2001). For a given 
galaxy of flux density Si 5 located at a given redshift z, 
we have computed S'15 for each one of the 105 template 
SED if it were located at the same redshift z. The tem- 
plate for which the computed 5' 15 was the closest to the 
observed S'15 was then used to compute the rest-frame 
15 /im luminosity, vL u [15 /jm] (open circles in Fig. 0b), 
after a normalization by a factor Sis/S'is. 

The rest-frame 15 /im luminosities estimated using the 
three techniques for the ISOCAM-HDFN galaxies are con- 
sistent within 20% (l-cr, see Fig.|b). Above z ~ 0.6, the 
first technique, which makes use of the SED of M 82, sys- 
tematically underestimates j/L„[15/Ltm], as compared to 
the other two techniques. This was to be expected since 
M 82 is a moderately luminous IR galaxy with a lower 
FIR over 6.75 /im luminosity ratio than more luminous 
galaxies (see Fig. ||d). 

In the following, we will use the multi-template tech- 
nique to estimate the IR luminosity and SFR of galaxies, 
because it is consistent with the correlations between MIR 
and FIR luminosities and provides a continuous interpo- 
lation between the 6.75, 12,15/xm and FIR luminosities 
based on observed spectra. Since the templates fit the cor- 
relation between Lir and ^L t/ [15/im], it is equivalent to 
compute vL u [lh /im] first and then use Eq. ( [l3|) which 
links Lir and vL u [lb /im] or to estimate directly Lir with 
the best template SED. 
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Fig. 8. Ltr [8- 1000 /im] and star formation rate (SFR) 
versus redshift and A rf for the ISOCAM-HDFN galaxies 
(S15 > 0.1 mJy). The 15 /im completeness limit of 0.1 
mJy is materialized with a dashed line. Only the left axi s 
is meaningful for the five AGNs (filled dots, Sect. 5.2). 
HDF_PS3_6a (no spectroscopic redshift) is not included. 



5.4- Luminosities and star formation rates 

In this section, we use the multi-template technique to di- 
rectly estimate the L m and SFR of the ISOCAM-HDFN 
galaxies. The results are shown in the Fig. ^ and Tabic || 
The l-cr error bars in Fig. ^| were computed from the dis- 
persion measured in the correlations of the Fig. |^. For 
example, for a galaxy located at z ~ 1, since A r f is close 
to 7 /Ltm, the error bar on Lir is the quadratic sum of the 
error bar on Eq. ( ]ll|) , which gives L15 from L7, and of the 
error bar on Eq. (|13|), which gives Lir from L15. 

The dashed line on Fig. || corresponds to the sensitivity 
limit of Si5= 0.1 mJy converted into a Lir(z) as in the 
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Fig. 9. a) 15 /im versus radio continuum (1.4 GHz) luminosities. Small filled dots: sample of 109 local galaxies from 
ISOCAM and NVSS. Filled dots with error bars: 17 HDFN galaxies in common between ISOCAM and VLA or WSRT 
(see text). Open dots with error bars: 7 CFRS-14 galaxies in common between ISOCAM and VLA (Flores et al. 1999). 
b) "q" parameter as a function of redshift for the HDFN and CFRS-14 galaxies. Plain line: median value of "q" 
(2.31q5). Dashed line: measured value from local galaxies (q— 2.34±0.01; Yun, Reddy & Condon 2001). Dotted line: 
l-cr error bar on "q" for the 24 galaxies. 



Fig.|. Table | summarizes the properties of the ISOCAM 
galaxies in each of these four galaxy categories plus the 
AGN category. 

About 75 % of the galaxies dominated by star forma- 
tion are either LIGs or ULIGs, while the remaining ~ 
25% are nearly equally distributed among either "star- 
bursts" (10 11 > (L IR /L Q ) > 10 10 ) or "normal" galaxies 
(Lir < 10 10 Lq). The contribution of each galaxy type to 
the 15 urn IGL (IGL\$) is given in the last column of the 
Table |3j. Luminous IR galaxies are responsible for about 
60% of IGL 15 while AGNs contribute to about 20%, the 
remaining 20 % being due to the normal and starburst 
galaxy types. This repartition will be used in Sect. || to 
compute the 140 /im IGL (IGL 1A q) from ISOCAM galax- 
ies. 

The mean IR luminosity of the ISOCAM-HDFN galax- 
ies is Lir ~ 6xl0 n Lq. It corresponds to a SFR ~ 100 
Mq yr _1 (from Eq. (0)). In the following section, we eval- 
uate the robustness of this result with an independent 
estimator of the IR luminosity: the radio luminosity. 



5.5. The radio-MIR correlation 

The radio continuum is also a tracer of the bolometric 
IR luminosity of star forming galaxies as a result of the 
tight correlation between IR and 1.4 GHz radio continuum 
in local galaxies (Condon 1992, Yun, Reddy & Condon 
2001; see Sect. [2~4| ). The origin of this correlation remains 
unclear, but it is generally assumed that massive stars are 
both responsible for the UV photons that heat dust before 
it radiates in the IR and for the synchrotron acceleration 



of electrons, producing the radio continuum, when they 
explode as supernovae. 

In Fig. |^a, we have plotted the 1.4 GHz and 15 /im lu- 
minosities of the 109 local galaxies (small filled dots) from 
the ISOCAM sample described in Sect. 4.2.1, which were 



detected in the NRAO VLA Sky Survey (NVSS, Condon 
et al. 1998). As expected both luminosities are correlated 
with each other since both the 1.4 GHz and 15/xm lumi- 
nosities are correlated with Lir. 

Half of the 35 star forming HDFN-ISOCAM galaxies 
with spectroscopic redshift (AGNs excluded) present flux 
densities larger than 40/xJy at 1.4 GHz (5-cr) in the VLA 
and WSRT catalogs from Richards (2000) and Garrett et 
al. (2000) respectively. Their rest-frame 15 /im luminosi- 
ties were computed in the Sect. 5.3. Their rest-frame 1.4 



GHz lu mino sities were computed assuming a power-law as 
in Sect. 2A: S u oc v~ a , where a= 0.8±0.15 as suggested 
for star forming galaxies in Yun, Reddy & Condon (2001). 
The rest-frame 15 /im and 1.4 GHz luminosities of these 
17 HDFN-ISOCAM galaxies are plotted as filled dots with 
error bars in the Fig. ||a. The error bars on vL v [\A GHz] 
were computed from the quadratic sum of the error bars 
on the measurement of the radio flux densities plus the 
error bar on a. Nine galaxies common in both catalogs 
from Richards (2000) and Garrett et al. (2000) present up 
to 40 % differences in their radio flux densities. For these 
galaxies, we used the mean value and included the differ- 
ence between both measurements in the error bars. 

We have also included 7 galaxies from the CFRS-14 
field (Flores et al. 1999), which rest-frame luminosities 
were computed with the same techniques (open dots). 
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The MIR and radio luminosities of this sample of 24 
distant dusty galaxies (z ~ 0.8) are also strongly corre- 
lated with each other. This is the hrst time that such cor- 
relation is found at these rcdshifts. The plain line in the 
Fig. ||a is a power-law fit to this correlation. For compar- 
ison, we have plotted in dashed line the correlation that 
one would expect from the combination of Eq. (|l^) (Ljr 
from L15) and the radio "q" parameter (see Sect. |2.4| ). 
The slope of the correlation observed for distant galax- 
ies is marginally steeper than this dashed line. It is also 
slightly steeper than the correlation found for local galax- 
ies (small filled dots), but the difference is too marginal 
with respect to the number of galaxies at high luminosities 
for detailed interpretation. 

In Fig. |b, the "<f parameter for the 24 ISOCAM- 
HDFN and CFRS-14 galaxies is plotted as a function of 
redshift. The "g" parameter is the logarithm of the ratio 
of Z/fir over L(1.4 GHz) (see Eq. (^)). Lir was estimated 
as described in Sect. 5.4, then converted into Lfir from 
Eq. (|J) . The median value of the "g" parameter for the 24 
galaxies is: q= 2.3j^g (plain and dotted lines in Fig. ||b), 
in perfect agreement with the local value of q= 2.34±0.01 
(Yun, Reddy & Condon 2001, dashed line in Fig. |b). 

This study shows that the IR luminosities estimated 
from the ISOCAM 15 /im flux densities using the MIR-FIR 
correlations are perfectly consistent with those estimated 
from the radio. Although it is not clear whether the radio 
versus IR correlation also applyies up to z ~ 1 , this result 
independently validates our estimate of the bolometric IR 
luminosity of the ISOCAM galaxies. 

5. 6. Evolution of the infrared volume emissivity with 
redshift 

We have seen that the excess of faint MIR sources in num- 
ber counts was due to the presence of distant (z ~ 1) lu- 
minous IR galaxies. As a consequence, the amount of star 
formation per comoving volume hidden by dust must have 
rapidly decreased from z~l to 0. Indeed, the ISOCAM- 
HDFN galaxies with 0.6 < z < 1.3 and L m > 10 11 Lq 
(LIGs and ULIGs; AGNs excluded) produce a 15 /im lu- 
minosity density of £15 [z ~ 1) ~ (60 ± 12) x 10 6 Lq 
Mpc~ 3 , while in the local universe luminous IR galax- 
ies only make C 15 {z ~ 0) ~ (11 ± 3) x 10 5 Lq Mpc~ 3 
(see Table |). The comoving luminosity density produced 
by luminous IR galaxies at 15 /im was therefore about 55 
times larger at z ~ 1 than in the local universe. The 15 /im 
luminosity densities at z ~ were calculated using the lo- 
cal luminosity function (LLF) from Xu et al. (1998). In 
order to estimate the contribution of luminous IR galax- 
ies, we used the Eq.([l3]) to convert the IR luminosity cut- 
off of Lir = 10 11 Lq into a 15 /im luminosity cut-off 
of L 15 ~ 10 10 Lq. The 15 /im luminosity density that 
we have computed is consistent with the one measured at 
12 /tm from IRAS, C z ^° ~ (1.7 ± 0.5) x 10 6 Lq Mpc" 3 , 
converted to 15 /mi using Eq.(n2h. 





15 /im 


8-1000 /tm 


0.44 /im 




All LIGs 


All LIGs 


All 


z ~0 


7.5±1.6 1.1±0.3 


120±30 8.8±1.6 


82±2 


z ~1 


61±12 


632±291 


240±100 



Table 4. Comoving luminosity density in units of 
10 6 Lq Mpc~ 3 in the local universe and at z ~ 1 (0.6 
< z > 1.3). Columns "AH": contribution of galaxies of all 
luminosities are considered. Columns "LIGs" : only galax- 
ies with Lir > 10 11 Lq. AGNs were excluded for the 
computation of the IR values. The B-band (0.44 /im) val- 
ues are from Lilly et al. (1996) at z ^0 and Connolly et 
al. (1997) at z ~ 1 (mean of the 0.5-1 and 1-1.5 redshift 
bins). 



If we now consider the bolometric IR luminosity (from 
8 to 1000 /im) of ISOCAM galaxies, as estimated in the 
Sect.|, we find that the comoving density of IR luminosity 
radiated by dusty starbursts was about (70 ± 35) times 
larger at z ~ 1 than today (computed from the LLF of 
Soifer et al. 1987). Since the IR luminosity is directly pro- 
portional to the extincted star formation rate of a galaxy, 
this means that the comoving density of star formation 
taking place in luminous IR galaxies was about (70 ± 35) 
times larger at z ~ 1 than today. In case of a pure density 
evolution proportional to (l+z) n , this would translate into 
a value of n ~ 6. However, we want to emphasize that we 
are only considering here the galaxies detected by ISO- 
CAM, not the full luminosity function. For comparison, 
the £?-band (0.44 /im) luminosity density was only about 
three times larger at z ~1 than today. 

The redshift evolution of the comoving density of IR 
luminosity is compared to the UV (2800A) one in Fig. 0a. 
Both wavelengths exhibit similar luminosity densities at 
low redshift but the IR rises faster than the UV and 
reaches a larger value at z ~ 1, implying that a much 
larger fraction of star formation was hidden by dust at 
z ~ 1 than today. 

Finally, we note that the projected density of galaxies 
detected in the £>-band in the HDFN (529 galaxies/'" with 
Bab < 29, Pozzetti et al. 1998) is 330 times greater than 
the projected density of ISOCAM galaxies (1.6 sources/' , 
with S15 > 0.1 mJy or AB(15 fim)< 18.9). However the 
ISOCAM galaxies produce a 15 /im IGL which is only 
twice lower than the B-band IGL (IGL B ~ 4.57tolr 
nW m~ 2 sr _1 , Madau & Pozzetti 2000). This confirms 
that the very luminous galaxies detected in the MIR ra- 
diate mostly in the IR. 



6. Contribution of ISOCAM galaxies to the 
cosmic infrared background 



Fig. |10| b summarizes the current estimates of the cos- 
mic background for A = 0.1 to 1500 /im. Lower limits to 
the EBL (filled dots) come from integrating galaxy fluxes 
(IGL) , while direct measurements of the EBL come from 
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Fig. 10. a) FIR (filled circles), MIR- 15 /mi (open circles) and UV-2800A luminosity density (in Lq Mpc~ 3 ) as a 
function of redshift. UV-2800A: open squares from Cowie, Songaila & Barger (1999); filled squares from Lilly et al. 
(1996). b) Integrated Galaxy Light (IGL, filled dots) and Extragalactic Background Light (EBL, open squares, grey 
area) from the UV to sub- millimeter (see Table || for the origin of the data). EBL measurements from COBE: 200- 
1500 urn EBL from COBE-FIRAS (grey area, Lagache et al. 1999), 1.25, 2.2, 3.5, 100, 140 /mi EBL from COBE-DIRBE 
(open squares, see Table. ||). IGL in the U,B,V,I,J,H,K bands from Madau & Pozzetti (2000). The upper end of the 
arrows indicate the revised values suggested by Bernstein et al. (2001, factor two higher). Our estimate of the 15 /tm 
IGL (2.4 ± 0.5 nW m~ 2 sr _1 ) is marked with a surrounded star. 6.75 /tm (ISOCAM-LW2 filter) IGL from Altieri et 
al. (1999, filled dot). Hatched upper limit from Mkn 501 (Stanev & Franceschini 1998). The arrow from 15 to 140 /tm 
indicates our computation of the 140 /tm IGL due to ISOCAM galaxies. 



the COBE experiments DIRBE (open squares) and FI- 
RAS (grey area, A > 200 /tm, Lagache et al. 1999). Table | 
gives the origin of the data used in Fig.[lO|b. Three inde- 
pendent studies obtain the same FIR to sub-millimeter 
EBL above A ~ 140 /tm. We did not consider here the 
EBL at A = 60 /tm (Finkbeiner, Davis & Schlegel 2000) 
because it strongly depends on the remains controversial 
(see Puget & Lagache 2001). 

Below A ~ 5 /tm, the cosmic background originates 
from direct stellar light, while above this wavelength, it 
comes from either stellar or AGN light reprocessed by 
dust. In the dust domain, i.e. from A ~ 5 to 1500 /tm, 
the cosmic background peaks around A ~ 140 /tm. It is 
not directly measured over this whole wavelength range, 
but below A ~ 40 /tm and down to about 2 /mi, an up- 
per limit is set by the observations of the TeV outburst of 
Mkn 501 (see Sect. §). 

The 15 /tm IGL (IGL 15 = (2.4±0.5) nW m~ 2 sr" 1 , see 
Sect. ||) is marked with a surrounded star on the Fig. |To|b. 
It is about ten times lower than the peak of the CIRB at 
Amax ~ 140 /tm. At the median redshift of the ISOCAM 
galaxies, z ~ 0.8 (Sect. 5.1), the rest-frame wavelengths 



corresponding to the observed 12-18 /tm (ISOCAM-LW3 
bandpass) and 140 /tm wavelengths, are 6.5-9.8 /tm and 
~ 80 /tm. These rest-frame wavelengths are close to the 
ISOCAM-LW2 (5-8.5 /tm) and IRAS (60, 100 /tm) bands. 
We have seen in the Fig. ||d that the luminosities of local 
galaxies in both bands correlate with each other. If this 



correlation remains valid up to z ~ 1, then it implies that 
we can compute the 140 /tm IGL due to ISOCAM galaxies. 
There is no obvious reason why the SEDs of galaxies at 
z ~ 1 would exhibit very different IR SEDs. For example, 
the SED of the extremely red object HR 10 [z ~ 1.44) is 
very similar to the one of Arp 220 normalized by a factor 
of - 4 (Elbaz et al. 2002). 



We propose here to use the multi-template technique, 
as in Sect. |], to compute the contribution of ISOCAM 
galaxies to the 140 /tm background. We will first sep- 
arate the relative contributions to IGL15 from ULIGs, 
LIGs, SBs, normal galaxies and AGNs as estimated in 
the Sect. |^. For each galaxy type, a median redshift, z me d, 
and the ratio TZ = vS v [l4Q/(l + z med )]/vS v [lb/(l + z med )] 
are estimated based on the HDFN sample (see Table 
Although a larger sample would improve the statistical 
reliability of this computation, this choice is justified by 
the fact that the bulk of the 15 /tm IGL is produced by 
galaxies with flux densities in the range Sis — 0.35-1 mJy 
as the HDFN and CFRS-14 galaxies, which share similar 
redshifts distributions. It is important to separate the con- 
tribution from galaxies of different luminosities since the 
IR versus 6.75 /tm luminosity correlation presents a slope 
larger than one (Fig. ||d). The contribution from AGNs is 
computed assuming that they all share the SED of NGC 
1068. This is a conservative hypothesis since NGC 1068 
exhibits the lowest FIR over MIR ratio that we know. 
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Table 5. Integrated Galaxy Light (IGL, upper half of the 
table) and Extragalactic Background Light (EBL, lower 
half of the table), vl u: in units of (*) nW m~ 2 sr _1 . (a) 
1-cr error bars. (b),(c) include the systematic error quoted 
by the authors. Ref: (1) Madau & Pozzetti (2000), (2) 
Altieri et al. (1999), (3) Matsuhara et al. (2000), (4) Puget 
et al. (1999), (5) Barger, Cowie & Sanders (1999), (6) 
Blain et al. (1999b), lensed galaxies Q), (7) Wright (2000), 
(8) Cambresy et al. (2001), (9) Dwek & Arendt (1998), 
(10) Gorjian, Wright & Chary (2000), see also Wright & 
Reese (2000), Wright (2000), (11) Finkbeiner, Davis & 
Schlegel (2000) and Hauser et al. (1998) for the 140 and 
240 (im values, (12) Lagache et al. (2000), (13) Lagache et 
al. (1999), Fixsen et al. (1998) 



The results are summarized in the Table We first 
separate the contribution of each type of galaxy to the 
15 £im IGL (Col.(3),(4)). The 140 /im IGL is then com- 
puted by converting each contribution to the 15 /im IGL 
into a contribution at 140 /im, by multiplying by TZ each 
time. We find a total of IGL li0 ~ (16±5) nW m~ 2 sr" 1 
(see Fig. |lO|b) , which corresponds to (65±35) % of the 
observed value from COBE-DIRBE (EBL U0 = (25±7) 
nW m" 2 sr" 1 , Hauser et al. 1998, Lagache et al. 2000, 
Finkbeiner et al. 2000). 



Table 6. Contribution of the ISOCAM galaxies to the 
140 /im extragalactic background light. Col.(l) galaxy 
type: AGNs and star forming galaxies above and below 
log{L m /L Q ) = 10, 11 and 12 (see Sect. z, in Col. (2), 
is the median redshift for each galaxy type. "%15", in 
Col. (3), is the fractional contribution of each galaxy type 
to IGLu (EqM)). Col.(4), IGL 15 for each galaxy type 
in nW nr 2 sr- V . Col. (5), TZ = vS v [U0 fJ,m]/vS v [15 fim}. 
For AGNs, we used the SED of NGC 1068. IGLuo, in 
Col. (4), is the 140 /im IGL produced by each galaxy type 
in nW m~ 2 sr _1 and %140, in Col. (5), is the fractional 
contribution to the 140 /im IGL for each galaxy type. 



Type 


z 


%15 


IGLxs 


K 


IGL140 


%140 


ULIGs 


1.2 


17 


0.4±0.1 


14±6 


5.7±2.8 


36 


LIGs 


0.8 


44 


1.1±0.2 


7±3 


7.7±3.6 


48 


SBs 


0.5 


9 


0.2±0.1 


4±1 


1.0±0.4 


6 


normal 


0.1 


11 


0.3±0.1 


4±1 


0.9±0.5 


6 


AGNs 


1.0 


19 


0.5±0.1 


1.4 


0.7±0.1 


4 


Total 


0.8 


100 


2.4±0.5 


7±3 


16.0±4.6 


100 



Hence the galaxies detected in ISOCAM deep surveys 
are found to be responsible for the bulk of the CIRB. 
About half of this 140 //m IGL is due to LIGs and about 
one third from ULIGs. The contribution of AGNs is esti- 
mated to be as low as ~ 4% of the 140 /im IGL. It could 
be as high as ~ 20 % if instead of NGC 1068, we had used 
the SED of NGC 6240. In this case, ISOCAM galaxies 
would produce a 1GL\^ ~ (19±5) nW m~ 2 sr _1 . If we 
had included the contribution from galaxies with flux den- 
sities between 30 and 50 /iJy only detected through lensing 
magnification (see Sect. |^), IGLi^q would be larger by an 
extra ~ 4 nW m~ 2 sr _1 . 

At the mean redshift of z ~ 0.6 where the bulk of 
the B-band IGL is produced (Pozzetti & Madau 2001), 
light was radiated in the UV in order to be observed at 
0.44 /im. As we have seen, the bulk of the 140 /im is also 
produced in the same redshift range. Hence the 140 //m 
and i?-band probe the UV emission from stars with and 
without dust extinction. The ratio of IGL 140 over IGLb 
should therefore give a rough approximation of the ratio 
of extincted over non extincted star formation around z ~ 
0.5-1. This ratio is close to 5, which suggests again that 
the bulk of the UV photons radiated by young stars in this 
redshift range was strongly affected by dust extinction. 
This value is consistent with the one obtained by Chary 
& Elbaz (2001) or Franceschini et al. (2001). 

7. Conclusions 

The cosmic IR background is a fossil record of the light 
radiated by galaxies since their formation. It peaks around 



A,. 



140 /im whereas the spectral energy distribution of 
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galaxies peaks above 60 /xm. This suggests that the bulk of 
the cosmic IR background is due to galaxies located below 
z ~ 1.3. We have shown that the best technique currently 
available to unveil dusty galaxies up to 2 ~ 1.3 is provided 
by the ISOCAM MIR extragalactic surveys. 

We have computed the contribution of ISOCAM galax- 
ies to the 15 fiin background, the 15 /im integrated galaxy 
light, and found a value of IGL 15 ~ (2.4±0.5) nW m~ 2 
sr _1 . This is about ten times below the cosmic background 
measured by COBE at A max ~ 140 (im: EBL 140 ~ (25±7) 
nW m~ 2 sr _1 . 

We have demonstrated that the MIR luminosities at 
6.75, 12 and 15 /im were correlated with each other and 
with the bolometric IR luminosity for local galaxies. This 
suggests that the contribution of ISOCAM galaxies to the 
CIRB can be computed from JGL15, unless distant galax- 
ies SEDs strongly differ from local ones. 

The redshift distribution of ISOCAM galaxies was 
measured from the spectroscopically complete sample of 
galaxies in the region of the HDFN. This redshift distribu- 
tion is consistent with the twice larger sample of ISOCAM 
galaxies detected in the CFRS fields CFRS-14 and CFRS- 
03 (Flores et al. 1999, 2002). At the median redshift of z ~ 
0.8, the observed 15 and 140 /jm wavelengths correspond 
to about 7mhi (ISOCAM-LW2 filter) and 80 fim (IRAS 
bands) in the rest-frame of the galaxies. Luminosities at 
both wavelengths are correlated (see Fig. ||d). If the cor- 
relations between MIR and FIR luminosities remain valid 
up to z ~ 1, then they can be used to compute IGL140. 
We have checked with a sample of galaxies detected both 
in the MIR with ISOCAM and in the radio with the VLA 
and WSRT, that the MIR- FIR and radio- FIR correlations 
are consistent up to z ~ 1. This comparison independently 
validates our estimate of the bolometric IR luminosity of 
the ISOCAM galaxies, although it is not clear whether the 
radio-FIR correlation works also up to z ~ 1. 

The fraction of active nuclei responsible for the 15 /im 
luminosity of ISOCAM galaxies was estimated from the 
deepest soft and hard X-ray surveys available at present 
by the XMM-Newton and Chandra X-ray observatories in 
the Lockman Hole and HDFN respectively (Fadda et al. 
2001). It was found that about (12±5) % of the ISOCAM 
galaxies are powered by an AGN and that the AGN con- 
tribution to IGL 15 was about (17±6) %. The AGN con- 
tribution to IGL140 was found to be as low as ~ 4 % as- 
suming that they all share the SED of the local Seyfert 
2, NGC 1068. This is a conservative choice since NGC 
1068 presents the flattest IR SED that we know. However, 
we note that the cosmic X-ray background (CXB) peaks 
around 30 keV (see Fig.l of Wilman, Fabian & Nulsen 
2000), while both XMM-Newton and Chandra were lim- 
ited to energies below 10 keV. It is therefore possible that 
a population of hard X-ray AGNs was missed by these sur- 
veys. But as claimed by the authors of these deep X-ray 
surveys, the bulk of the CXB had been resolved into indi- 
vidual galaxies in the Lockman Hole and HDFN images. 



Moreover, using estimates of the present comoving density 
of black holes, Madau & Pozzetti (2000) calculated that 
less than 20 % of the CIRB could be due to dusty AGNs. 

For the remaining star forming galaxies, we used a li- 
brary of template SEDs, reproducing the MIR-FIR cor- 
relations, to compute their IR luminosity and contribu- 
tion to IGL14Q. We find that LIGs and ULIGs produce 
about 60% of IGL15. The comoving density of IR lumi- 
nosity produced by these luminous IR galaxies was about 
(70±35) times larger at z ~ 1 than today, while in the 
same redshift interval the i?-band or UV luminosity den- 
sities only decreased by a factor ~ 3. Since the IR lumi- 
nosity measures the dusty star formation rate of a galaxy, 
this also implies that the comoving density of star forma- 
tion, due to luminous IR galaxies, decreased by a similar 
factor in this redshift range, i.e. much more than expected 
by studies at other wavelengths affected by dust extinc- 
tion. This indicates that a large fraction of present day 
stars were formed during a dusty starburst event. 

We estimate a contribution of ISOCAM galaxies to the 
peak of the CIRB at X max ~ 140 ^m of (16±5) nW mT 2 
sr _1 as compared to the measured value of (25±7) nW 
m~ 2 sr _1 from COBE. This study therefore suggests that 
the ISOCAM galaxies are responsible for the bulk of the 
CIRB. 

We have started a systematic spectroscopic follow-up 
of these galaxies with the aim of studying their physi- 
cal properties and the origin of their large SFR. Frances- 
chini et al. (2001) estimated their baryonic masses to be of 
the order of (M) ~ 10 11 Mq by fitting their optical and 
near-IR luminosities with template SEDs (from Silva et al. 
1998) and assuming a Salpeter initial mass function (from 
0.15 to IOQMq). Their colors are similar to field galax- 
ies of similar magnitudes (Cohen 2001), hence they could 
not have been selected on the basis of their optical colors. 
The technique which consists in using the spectral slope 
in the UV domain to correct luminosities from extinction 
(Meurer, Heckman & Calzetti 1999) fails to differentiate 
the luminous dusty galaxies detected with ISOCAM from 
other field galaxies in the HDFN+FF (Cohen 2001). This 
was to be expected since this technique only works for 
galaxies with L m < 4 x 10 11 Lq (Meurer et al. 2000) 
while most ISOCAM galaxies are more luminous than this 
threshold. A property of the ISOCAM galaxies that may 
give a hint on their origin is their association with small 
groups of galaxies. A preliminary study of 22 ISOCAM- 
HDFN galaxies by Cohen et al. (2000) found that nearly 
all ISOCAM galaxies belonged to small groups, while the 
fraction of field galaxies with similar optical magnitudes 
belonging to such groups was 68 %. The study of the full 
sample of ISOCAM-HDFN galaxies by Aussel et al. (2001) 
shows that all of them belong to physical groups, hence 
suggests that dynamical effects such as merging and tidal 
interactions are responsible for their large SFR. 

The Space IR Telescope Facility (SIRTF) will soon 
provide a powerful insight on the FIR emission of distant 
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dusty starbursts and its 24 /im band, less affected by con- 
fusion, should prolongate their detection up to z ~ 2.5. 
Extending the redshift range surveyed by ISOCAM to 
such redshifts is particularly important to measure the di- 
rect IR emission of the distant population of Lyman break 
galaxies, whose dust extinction remains highly uncertain 
(Steidel et al. 1999). However, the direct measurement of 
the MIR and FIR emission of distant galaxies will only 
become possible with the launch of the Next Generation 
Space Telescope (NGST) and Herschel (FIRST) satellite 
scheduled for 2009 and 2007. 

Acknowledgements. DE wishes to thank the American Astro- 
nomical Society for its support through the Chretien Interna- 
tional Research Grant and Joel Primack and David Koo for 
supporting his research through NASA grants NAG5-8218 and 
NAG5-3507. We wish to thank H.Flores and A.Blain for fruitful 
comments and A.Boselli for sharing his ISOCAM catalog with 
us. This research has made use of the NASA/IPAC Extragalac- 
tic Database (NED) which is operated by the Jet Propulsion 
Laboratory, California Institute of Technology, under contract 
with the National Aeronautics and Space Administration. 

References 

Allamandola, L.J., Tielens, A.G.G.M., Barker, J.R. 1989, ApJS 
71, 733 

Almaini, O., Lawrence, A., Boyle, B.J. 1999, MNRAS 305, L59 
Altieri, B., Metcalfe, L., Kneib, J. P., et al. 1999, A&A 343L, 
65 

Anantharamaiah, K.R., Viallefond, F., Mohan, NR., Goss, 

W.M., Zhao, J.H. 2000, ApJ 537, 613 
Andriesse, CD., 1978, A&A 66, 169 

Aussel, H., Cesarsky, C.J., Elbaz, D., Starck, J.L. 1999, A&A 
342, 313 

Aussel, H., Coia, D., Mazzei, P., De Zotti, G, Franceschini, A. 

2000, A&AS 141, 257 

Aussel, H., et al. 2001, in preparation 

Barger, A. J., Cowie, L.L., Sanders, D.B. 1999, ApJ 518, L5 
Bernstein, R. A., Freedman, W . L., Madore, B.F. 2001, ApJ, 

accepted ( ^stro-ph/011219lj ) 
Blain, A.W., Smail, I., Ivison, R.J., Kneib, J.-P. 1999a, MN- 
RAS 302, 632 

Blain, A.W., Kneib, J.-P, Ivison, R.J., Smail, I. 1999b, ApJ 
512, L87 

Boselli, A., Lequeux, J., Sauvage, M., et al. 1998, A&A 335, 53 
Brandt, W.N., Hornschemeier, A.E., Alexander, D.M., et al. 

2001, AJ 122, 1 

Brandt, W.N., Alexan der, P.M., Hornsc hemeier, A.E., et al. 

2002, AJ, in press (|astro-ph/0108404|) 



Cambresy, L., Reach, W.T., Beichman, C.A., Jarrett, T.H. 

2001, ApJ 555, 563 
Cesarsky, D., Lequeux, J., Abergel, A., et al. 1996a, A&A 315, 

L305 

Cesarsky, C.J., Abergel, A., Agnese, P., et al. 1996b, A&A 315, 
L32 

Charmandaris, V., Laurent, O., Mirabel, I.F., et al. 1999, 

Ap&SS 266, 99 
Chary, R.R., Elbaz, D. 2001, ApJ 556, 562 
Cohen, J.G, Hogg, D.W., Blandford, R., et al. 2000, ApJ 538, 

29 



Cohen, J.G. 2001, AJ 121, 2895 

Coleman, G. D., Wu, C.-G, Weedman, D. W. 1980, ApJS 43, 
393 

Condon, J.J., Anderson, M.L., Helou, G. 1991, ApJ 376, 95 
Condon, J.J., 1992, ARA&A 30, 575 

Condon, J. J., Cotton, W.D., Greisen, E.W., et al. 1998, AJ 
115, 1693 

Connolly, A.J., Szalay, A.S., Dickinson, M.E., et al. 1997, ApJ 
486, Lll 

Cowie, L.L., Songaila, A., Barger, A.J. 1999, ApJ 118, 603 
Dale, D.A., Silbermann, N.A., Helou, G, et al. 2001, AJ 120, 
583 

Desert, F.-X., Boulanger, F., Puget, J.-L. 1990, A&A 237, 215 
de Vaucouleurs, G, de Vaucouleurs, A., Corwin, H.G.Jr., et al. 

1991, Third reference catalogue of Bright galaxies, (New 

York : Springer- Verlag) (RC3) 
Dole, R, Gispert, R., Lagache, G, et al. 2001, A&A 372, 364 
Dwek, E., Arendt, R.G, Hauser, M., et al. 1998, ApJ 508, 106 
Dwek, E., Arendt, R.G. 1998, ApJ 508, L9 
Eales, S., Lilly, S., Webb, T., et al. 2000, AJ 120, 2244 
Elbaz, D., Cesarsky, C.J., Fadda, D., et al. 1999, A&A 351, 

L37 

Elbaz, D., Flores, R, Chanial, P., et al. 2002, A&A 381, LI 
Fadda, D., Elbaz, D., Due, P.-A., et al. 2000, A&A 361, 827 
Fad da, P., Flores, H, Hasinger, G, et al. 2001, A&A in press 
(|astro-ph/0111412j) 



Finkbeiner, P.P., Pavis, M., Schlegel, P.J. 2000, ApJ 524, 867 
Fischer, J., Luhman, M.L., Satyapal, S. 1999, Ap&SS 266, 91 
Fixsen, P. J., Pwek, E., Mather, J.C., Bennett, C.L., Shafer, 

R.A. 1998, ApJ 508, 123 
Flores, R, Hammer, F., Thuan, T.X., et al. 1999, ApJ 517, 

148 

Flores, R, et al. 2002, in preparation 

Forster Schreiber, N.M., Laurent, O., Sauvage, M., et al. 2000, 

A&A, in preparation 
Forster Schreiber, N.M., Genzel, R., Lutz, P., Kunze, P., 

Sternberg, A. 2001, ApJ 552, 544 
Franceschini, A., Aussel, R, Cesarsky, C.J., Elbaz, P., Fadda, 

P. 2001, A&A 378, 1 
Franceschini, A., Fadda, P., Cesarsky, C.J., Elbaz, P., Flores , 

R, Granato, GL. 2002, ApJ accepted (kstro-ph/0111413j) 



Garrett, M.A., de Bruyn, A.G, Giroletti, M., et al. 2000, A&A 
361, L41 

Gorjian, V., Wright, E. L., Chary, R.R. 2000, ApJ 536, 550 
Genzel, R., Lutz, P., Sturm, E., et al. 1998, ApJ 498, 579 
Genzel, R., Cesarsky, C.J. 2000, ARA&A 38, 761 
Gispert, R., Lagache, G, Puget, J.L. 2000, A&A 360, 1 
Haas, M., Klaas, U., Midler, S.A.H., Chini, R., Coulson, I. 

2001, A&A 367, L9 
Hasinger, G, Altieri, B., Arnaud, M., et al. 2001, A&A 365, 

L45 

Hauser, M.G, Arendt, R.G, Kelsall, T., et al. 1998, ApJ 508, 
25 

Helou, G, Khan, I.R., Malek, R., Boehmer, L. 1988, ApJS 68, 
151 

Hornschemeier, A.E., Brandt, W.N., Garmire, G.P., et al. 2001, 
ApJ 554, 742 

Hughes, P. R, Gear, W.K., Robson, E.I. 1994, MNRAS 270, 
641 

Hughes, P. R, Serjeant, S., Punlop, J., et al. 1998, Nature 
394, 241 



The bulk of the cosmic infrared background resolved by ISOCAM 



19 



Ikebe, Y., Leighly, K., Tanaka, Y., et al. 2000, MNRAS 316, 
433 

Jones, A. P., d'Hendecourt, L. 2000, A&A 355, 1191 
Kennicutt, R.C.Jr 1998, ARA&A 36, 189 

Kessler, M., Steinz, J., Anderegg, M., et al. 1996, A&A 315, 
L27 

Klaas, U., Haas M., Heinrichsen, I., Schulz, B. 1997, A&A 325, 
L21 

Lagache, G., Abergel, A., Boulanger, F., Desert, F.-X., Puget, 

J.-L. 1999, A&A 344, 322 
Lagache, G., Haffner, L.M., Reynolds, R.J., Tufte, S.L. 2000, 

A&A 354, L247 
Laurent, O., 1999, PhD thesis, Universite de Paris XI. 
Laurent, O., Mirabel, I.F., Charmandaris, V., et al. 2000, A&A 

359, 887 

Le Floc'h, E., Mirabel, I.F., Laurent, O., et al. 2001, A&A 367, 
L487 

Leger, A., Puget, J.L., 1984, A&A 137, L5 
Lemke, D., Klaas, U., Abolins, J., et al. 1996, A&A 315, L64 
Lilly, S.J., Lefevre, O., Hammer, F., Crampton, D. 1996, ApJ 
460, LI 

Madau, P., Pozzetti, L. 2000, MNRAS 312, L9 
Matsuhara, H., Kawara, K., Sato, Y., et al. 2000, A&A 361, 
407 

Metcalfe, L. 2000, In "The Extragalactic Infrared Background 
and its Cosmological Implications", IAU symposium 204, 
p.18 

Meurer, G.R., Heckman, T.M., Calzetti, D., 1999, ApJ 521, 64 
Meurer, G.R., Heckman, T.M., Seibert, M., et al. 2000, To ap- 
pear in Proc. JD9, IAU Manchester 2000, 'Cold Gas and 
Dust at High Redshift' ed. D.J.Wilner, H ighlights of As- 
tronomy Vol. 12. PASP (kstro-ph/001120l|) 



Mirabel, I.F., Vigroux, L., Charmandaris, V., et al. 1998, A&A 
333, LI 

Moutou, C, Verstraete, L., Leger, A., Sellgren, K., Schmidt, 
W. 2000, A&A 354, L17 

Okumura, K., 1998, ESA ISOCAM PSF Report, 
www.iso.vilspa.esa.es /users / expLlib/C AM / psf_rep.ps. gz 

Pozzetti, L., Madau, P., Zamorani, G., Ferguson, H.C., 
Bruzual, G.A. 1998, MNRAS 298, 1133 

Pozzetti, L., Madau, P. 2001, In "The Extragalactic In- 
frared Background and its Cosmological Implications", 
I AU Symposium 204, eds. M.Harwit and M.G.Hauser 
( tMtro-ph/0011359| ) 

Puget, J.-L., Leger, A., 1989, ARA&A 27, 37 

Puget, J.L., Abergel, A., Bernard, J. P., et al. 1996, A&A 308, 
L5 

Puget, J.L., Lagache, G., Clements, D., et al. 1999, A&A 345, 
29 

Puget, J.L., Lagache, G. 2001, In "The Extragalactic Infrared 
Background and its Cosmological Implications" , IA U Sym 



pos ium 204, eds. M.Harwit and M.G.Hauser, p. 19 (astro 
ph/|0101105|) 



Renault, C, Barrau, A., Lagache, G., Puget, J.-L. 2001, A&A 
371, 771 

Richards, E.A. 2000, ApJ 533, 611 

Rigopoulou, D., Spoon, H.W.W., Genzel, R., et al. 1999, AJ 
118, 2625 

Roussel, H., Vigroux, L., Bosnia, A., et al. 2001, A&A 369, 473 
Rush, B., Malkan, M.A., Spinoglio, L. 1993, ApJS 89, 1 
Sanders, D.B., Mirabel, I.F. 1996, ARA&A 34, 749 



Serjeant, S., Oliver, S., Rowan-Robinson, M., et al. 2000, MN- 
RAS 316, 768 

Silva, L., Granato, G.L., Bressan, A., Danese, L. 1998, ApJ 
509, 103 

Smail, I., Ivi son, R.J., Blain, A W, Kneib, J.-P. 2001, MNRAS 
in press (astro-ph/011210C ) 



Smith, C.H, Aitken, D.K, & Roche, P.F. 1989, MNRAS 241, 
425 

Soifer, B.T., Sanders, D.B., Madore, B.F., et al. 1987, ApJ 320, 
238 

Stanev, T., Franceschini, A. 1998, ApJ 494, L159 

Steidel, C.C., Adelberger, K.L., Giavalisco, M, Dickinson, M, 

Pettini, M. 1999, ApJ 519, 1 
Telesco, CM., Becklin, E.E, Wynn- Williams, C.G., Harper, 

D.A. 1984, ApJ 282, 427 
Tran, Q.D., Lutz, D., Genzel, R., et al. 2001, ApJ 552, 527 
Xu, C, Hacking, P.B., Fang, F., et al. 1998, ApJ 508, 576 
Xu, C. 2000, ApJ 541, 134 

Xu, C, Lonsdale, C.J., Shupe, D.L, O'Linger, J, Masci, F. 

2001, ApJ 562, 179 
Vigroux, L., Mirabel, F., Altieri, B., et al. 1996, A&A 315, L93 
Wilman, R.J., Fabian, A.C., Nulsen, P.E.J. 2000, MNRAS 319, 

583 

Wright, E.L. 2000, ApJ 553, 538 

Wright, E.L, Reese, E.D. 2000, ApJ 545, 43 

Yun, M.S., Reddy, N.A, Condon, J.J. 2001, ApJ 554, 803 



